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Abstract 

^ ■ To investigate the possibility of a ghost-antighost condensate the coupled Dyson-Schwinger equa- 
,— 1 ■ tions for the gluon and ghost propagators in Yang-Mills theories are derived in general covariant 
. gauges, including ghost-antighost symmetric gauges. The infrared behaviour of these two-point 
CO \ functions is studied in a bare-vertex truncation scheme which has proven to be successful in Lan- 
dau gauge. In all linear covariant gauges the same infrared behaviour as in Landau gauge is found: 
The gluon propagator is infrared suppressed whereas the ghost propagator is infrared enhanced. 
^ ■ This infrared singular behaviour provides indication against a ghost-antighost condensate. In the 
O . ghost-antighost symmetric gauges we find that the infrared behaviour of the gluon and ghost prop- 
. . , agators cannot be determined when replacing all dressed vertices by bare ones. The question of a 
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BRST invariant dimension two condensate remains to be further studied. 
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INTRODUCTION 



A large body of experimental data supports the general believe that Quantum Chromo- 
dynamics (QCD) is the correct theory of strong interactions. Nevertheless we are left with 
the task of understanding the physics of hadrons, and hereby in particular the mechanisms 
of confinement and spontaneous breaking of chiral symmetry. Gaining such an insight re- 
quires reliable non-perturbative treatments of QCD. Hereby Monte Carlo lattice calculations 
provide a rigorous non-perturbative approach to QCD. They have the advantage of fully re- 
specting gauge invariance independently of the size of the lattice used. On the other hand, 
the extraction of the continuum values of physical observables from the lattice data requires 
a careful study of the scaling regime. The observed scaling behaviour, however, will be in 
general contaminated by finite size effects. With respect to studies of the confinement mech- 
anisms this is problematic: As infrared singularities are expected to occur in QCD there is 
definite need for a continuum-based non-perturbative approach. 

To this end we note that the Schwinger-Dyson equations of QCD can address directly 
the infrared region. They provide genuine non-perturbative information and are at the same 
time fully formulated in the continuum theory. Such an approach is, however, less rigorous 
than lattice calculations in the sense that truncations of the tower of coupled equations are 
necessary in practical calculations. Justifications for such truncations can be given on the 
basis of general principles like e.g. a restriction to the first Gribov region, see ref. [I] and 
references therein. Nevertheless, the validity of the employed truncation is finally judged by 
comparing its results with either the results of Monte Carlo calculations or experiments. The 
latter is easily possible as the Schwinger-Dyson approach has been successfully applied to the 
description of hadron phenomenology, see e.g. the recent reviews refs. [E |3j] and references 
therein. Furthermore, despite recent progress by improved lattice algorithms, and despite 
the increasing computer time available for lattice calculations, including dynamical fermions 
is exceedingly cumbersome and finite baryon densities are hardly accessibly in realistic SU(3) 
lattice simulations. On the other hand, dynamical fermions and finite baryon densities can 
be relatively easily treated in the Schwinger-Dyson approach to QCD. 

In recent years the fundamental Schwinger-Dyson equations of SU(N) Yang-Mills theories 
have been solved explicitly in certain approximations yielding gluon and ghost propagators 
[E E3 IE 0, IE El- I n these calculations, carried out in Landau gauge, vertex functions 
constructed from appropriate Slavnov- Taylor identities as well as bare vertices have been 
employed. The results proved to b e q ualitati vely similar among each other and agree well 
with recent lattice calculations IE 11, IE IE 3 f° r both, the gluon and ghost propagator. 
The common, though gauge dependent, result of both approaches is an infrared suppressed 
gluon propagator and an infrared enhanced ghost propagator. Furthermore, the inclusion of 
dynamical quarks does not alter the infrared behaviour of gluon and ghost propagators and 
leads to only slight modifications for non- vanishing momenta for the number of light flavours 
Nf < 3 [15]. These results especially imply that the ghosts take the role of the long range 
correlations in the theory. Such a behaviour is in accordance with the Gribov-Zwanziger 
horizon condition, see e.g. ref. |7| and references therein, and the Kugo-Ojima confinement 
criterion, which in Landau gauge includes the statement that the ghost propagator should 
be more singular than a simple pole [la ]. 

The central assumption in the Kugo-Ojima confinement scenario is the invariance of 
the measure of the functional integral under BRS transformations and the existence of a 
nilpotent BRS-operator • The most general Lorentz invariant and globally gauge invariant 
Lagrangean of dimension four that can be constructed under this assumption has been 
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derived in ref. 18]. In addition to the structure appearing in ordinary linear covariant 
gauges, the Lagrangean contains a second gauge parameter which controls the symmetry 
of the Lagrangean under ghost- ant ighost interchange. Furthermore a four ghost interaction 
term is present. We will use this Lagrangean as the starting point of our investigation. 

Our main interest in this paper will be to explore the situation in these general covari- 
ant gauges. Away from the Landau gauge limit the connection between the Kugo-Ojima 
confinement criterion and the infrared behaviour of the ghost dressing function is far from 
obvious. In particular, the question might arise whether it is possible that the infrared 
dominant role of the ghost dressing function, seen in the Landau gauge, is assumed by other 
degrees of freedom like the longitudinal gluons in other covariant gauges. As a matter of fact, 
infrared dominance of longitudinal gluons is seen if stochastic quantization is used instead 
of the Faddeev-Popov quantization |l9j ]. Furthermore, calculations based on many-body 
techniques provide evidence that in Coulomb gauge (employing the usual Faddeev-Popov 
quantization) the ghosts and the Coulomb gluons are both infrared enhanced This 
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latter picture for Coulomb gauge QCD obtains (at least partial) support from lattice 
and renormalization group calculations 2^| . Care has, however, to be taken as the Coulomb 
gauge limit is highly non-trivial, see ref. e.g. j23[. On the other hand, the benefit of Coulomb 
gauge is obvious. The time-time component of the gluon propagator and the heavy quark 
potential fulfill a strictly valid inequality (22I 24j with the Coulomb string tension being 
several times larger than the asymptotic one Even more important, quark confinement 
directly results from infrared enhanced Coulomb gluons, see e.g. refs. |2g, |27[ and references 
therein. Instead of exploring the correlation functions in non-covariant gauges we will in this 
paper study Green's functions in covariant albeit non-linear gauges. 

Ghost-antighost symmetric gauges are of special interest when investigating the possibility 
of a BRST invariant condensate of dimension two in QCD. Such condensates occur in the 
operator product expansion of the gluon propagator [2^, 30|], bear some relation to the 
Gribov problem j3jj , may result in gluon mass generation |32l | and may be important for 
confinement in gener al I33L l34l . Hereby it has been clarified recently that these condensates 
are highly non-local [3jJ |36| and that they are only BRST invariant after eliminating the 
Nakanishi-Lautrup field via its equation of motion [37}. This kind of restricted BRST 
invariance has been called 'on-shell BRST invariance' and can be related to a residual gauge 
symmetry after gauge fixing. 

The solutions of the gluon and ghost Dyson-Schwinger equations in Landau gauge pro- 
vide a somewhat different picture: Whereas the operator product expansion of the gluon 
propagator requires such a dimension two condensate its interpretation with respect to a 
gluon mass is made impossible by the gluon propagator's infrared behaviour D(p 2 = 0) = 
instead of D(p 2 = 0) = 1/m 2 . Also the highly infrared singular ghost propagator excludes a 
ghost mass and/or a ghost-antighost condensate. Therefore the question arises whether in 
general ghost-antighost symmetric gauges the infrared behaviour of the propagators can be 
interpreted in terms of gluon and ghost "masses" . 

This paper is organized as follows: In sect. 2 we summarize some properties of the general 
Lagrangean given in ref. and outline the derivation of the coupled set of Dyson-Schwinger 
equations (DSEs) for the ghost and gluon propagators. As the Lagrangean contains a four- 
ghost interaction a rich structure in the ghost DSE emerges which closely resembles the one 
already present in the gluon equation of ordinary linear covariant gauges. In sect. 3 we 
employ a truncation scheme that has proven to be successful in Landau gauge and study 
in particular the infrared behaviour of the ghost and gluon dressing functions for general 
values of the two gauge parameters. Furthermore, we show that in the ghost-antighost 
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symmetric gauges the contributions of the genuine two-loop terms (generalized squint and 
sunset diagram) in the gluon and the ghost DSEs must be properly taken into account in 
the infrared. In the linear covariant gauges no such terms are present in the ghost DSE, and 
selfconsistent results can be obtained assuming the two-loop terms in the gluon equation to 
be subleading in the infrared • In general ghost-antighost symmetric gauges, on the other 
hand, the bare-vertex truncation is insufficient to clarify the infrared behaviour of the gluon 
and ghost propagators. In sect. 4 we will provide numerical solutions for the DSEs in the 
Landau gauge limit of the ghost-antighost symmetric case of the Lagrangean and recover 
the solutions found in from a different direction in two dimensional gauge parameter 
space. In the last section we give our conclusions. Technical details are deferred into four 
appendices. 



THE DYSON SCHWINGER EQUATION FOR THE GHOST PROPAGATOR 
Renormalized double BRS symmetry 

The most general Lagrangean of dimension four that is Lorentz invariant, globally gauge 
invariant, invariant under BRST- and anti-BRST-transformations, hermitean and omitting 
topological terms, is 



c = \f% + £40. + f (i-f)^(cxc)'- i\n^c -i(i-f) V-zv. (i) 

The field strength tensor and the covariant derivative are defined as 

F; v = d,A a u - d u Al - gt hc A\A c u 

Df = d,6 ab + gf abc A% (2) 

and the abbreviation (cx c) a = gf abc c b c c is used. Note that both ghost and antighost fields, 
c and c, resp., are chosen to be hermitean, c' = c and = c . This is necessary to maintain 
the hermiticity of the Lagrangian for all values of the gauge parameters A and a, see e.g. 
[l7l ] and references therein. Furthermore we work in Euclidean space-time. 

From the two gauge parameters of the Lagrangian the first one, A, is the usual parameter 
of linear covariant gauges, whereas the second one, a, controls the symmetry properties of 
the ghost content. For the cases a = and a = 2 one recovers the usual Faddeev-Popov 
Lagrangian and its mirror image, respectively, where the role of ghost and antighost have 
been interchanged. For the value a = 1 the Lagrangian is completely symmetric in the ghost 
and antighost fields. 

In ref. Q it has been shown that the S-matrix of the theory is invariant under variation 
of the gauge parameters A and a. Therefore gauge invariance of physical observables is 
ensured. One-loop calculations confirm in particular the independence of the first nontrivial 
coefficient of the (3 function from the gauge parameters. 

Furthermore, the existence of a renormalized BRS-algebra has been proven thus the 
theory given by (P) is multiplicatively renormalizable. From one-loop calculations one finds, 
that the Faddeev-Popov values of the gauge parameters, a = and a = 2, are fixed points 
under the renormalization procedure. The same is true for the ghost-antighost symmetric 
case a = 1. The case of Landau gauge, A = 0, corresponds to a fixed point as well, because 
the constraint d^A 11 = is not affected by a rescaling of the gluon field. 
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To be specific the renormalized BRS (s r ) and anti-BRS (s r ) transformations are given by 



i . ~ i 

-(cxc), s T c=-Z 1 - 



s r c = B - (cxc), s r c = -B - (1 - t;)Zi ( c x c ) > ^ 

*B = -fZx (c x B) S r S = -(1-1)^(5x5) 

"f (1 " f )^i 2 ((c x c) x c) , +f (1 - f )^ ((c x c) x c 

1 /2 

Here D r = d-Z^Z g (A x ) is the covariant derivative in the adjoint representation, with 
color and Lorentz indices suppressed. Note that the Nakanishi-Lautrup auxiliary field B 
can be eliminated from the BRS-transformations by using its equation of motion. The cor- 
responding BRS-transformations are called 'on-shell'. Note furthermore that the application 
of the BRS-operator s r (s r ) on a field increases (decreases) the ghost number by +1 (— 1), 
thus we can assign the value Npp = +1 (Npp = —1) to the (anti-)BRS-operator itself. The 
BRS-operator and the anti-BRS-operator are nilpotent and related by s r s r + s r s r = 0. These 
properties are, however, lost when considering 'on-shell' BRS-transformations. 

The Maurer-Cartan conditions, in addition to the forms of s r c and s r c, for ghosts c and 
anti-ghosts c in a ghost anti-ghost symmetric formulation thereby require 

s r c + s r c + Zx (c x c) = . (4) 

The correspondence between the bare Lagrangean and its renormalized version including 
counterterms is given by the following rescaling transformations 

Al -> y^A", c a c b -> Z 3 c a c b , B a -> B a /y% 
9 -> Z g g, a — ► Z a a, A -> Z X X, (5) 

where five independent renormalization constants Z3, Z3, Z g , Z a and Z\ have been intro- 
duced. Furthermore four additional renormalization constants are related to these via 
Slavnov-Taylor identities, 

Z\ = Z g Z\ , Zi = ZgZ^zJ , Z 4 = Z 2 g Z\, Z 4 = Z' l g Z\. (6) 

Note, however, that contrary to standard Faddeev-Popov gauges Z\ 7^ Z3, e.g., at one-loop 
(MS scheme), one has (HI 

" " ?21 :^(1-SA. (7) 



167r 2 e C 2 V 2 

The gauge fixing part of the Lagrangean (P) can be written in the following three equiv- 
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alent ways, 



Cqt = l;—^s r s r (z 3 AA + iZ x Z 3 aXcc) + ^(1 - a)^- s r (cs T c) 
2z 3 Z 3 V ) Z 3 y '2 



lS r C 



.Zx\ 
Z 3 2 

Z\ A 2 Z\ —r, a ,^ a, A,_ v , . .~ a s _„„ a 



( 9 ) 



= wdA +y 3 2 b2 + t z ' i (1 ~ r 2 ( " x c)2 + * Zs((1 ~ i } 59£V c + i 5 ^ • 

(10) 

This is verified by direct calculation via the transformations defined in Eqs. 0- In the form 
of Eq. (fTU|) the gauge fixing Lagrangean shows that the renormalization constants introduced 
in Q correspond to the replacements of bare by renormalized quantities as given above. 
We may rewrite the gauge fixing Lagrangean of Eq. (jl(Jj) once more, 

Lg T = iBdA + ^B 2 + §^ 2 £(1 - £(g x c) 2 + %Z 3 \ (cdD r c + cD r dc) 
Zj 3 Z Zj 3 A A l l 

+ iZ x {l -a) ~dA{cx c) . (11) 

This emphasizes the role of the gauge parameter a. In this form, the only term not symmetric 
under Faddeev-Popov conjugation, c — > c and c — > — c, is the last one (which is anti- 
symmetric w.r.t. Faddeev-Popov conjugation). It vanishes for a—1. With the current (real) 
hermiticity assignment for ghost and anti-ghost fields the Lagrangean is hermitean for all 
a, and it reduces to the standard Faddeev-Popov form for a = 0. We could also introduce 
hermitean adjoint ghost and anti-ghost fields, with the assignment = c, via the Caley map 
3|- This would then lead to 



C% F = iBdA + I^B 2 - ^-Zl |(1 - |) ^(c x c) 2 + Z 3 \{cdD r c + cD r dc) 

+ Z 1 (l-a)-(c{dAxc)-c{dAxc)) . (12) 

While this form of the Lagrangean, which we will not use further herein, is still hermitean 
it no longer reduces to the form of standard Faddeev-Popov theory for a = 0. Thus the 
Faddeev-Popov Lagrangean is only consistent with hermiticity for the choice of real ghost 
fields [13] . With complex conjugate ghost and anti-ghost fields, additional terms for a = 
survive (which are absent in standard Faddeev-Popov gauges). Only for a = 1 both versions, 
with hermitean real or complex conjugate ghost pairs, have the same Lagrangean and may 
be interchanged arbitrarily. 



Ghost and Anti-Ghost Dyson Schwinger equations 

Without invariance under Faddeev-Popov conjugation, i.e. without ghost anti-ghost sym- 
metry (a = 1 or A = 0), we have separate ghost and anti-ghost DSEs which are not identical. 
Consider the following representations of the ghost (anti-ghost) derivatives of the action (for 



6 



brevity we indicate by subscripts the space-time arguments of fields), 

|^ = Z 3 i(dD r c) a x + ^(l--)(cxB)» - ^Z 1 2 -(l--)((cx C )x C f 



= Z 8 i(0A.c)£ - = -<a r (5X; - i-£-B*) , (13) 

g = Z 8 i(^rc)2 + f (c x S) 2 + ^ Z^(l - |)((c x c) x c): 

= Z 3 i(dD r c) a x - ^s r B a x = -is r (dA% - i^Bl) . (14) 

^3 ^3 



The two DSEs then follow readily from 

<§4> = <<il§> = < i5 > 



Of course, they are related by Faddeev-Popov conjugation C FP which interchanges the two. 
In particular, 

C FP c = c , C FP c = — c , C F p£> = B + Zi(l — a)(c x c) , C FP y4 = A . (16) 

The transformation of the Nakanishi-Lautrup l?-field follows from compatibility with 
BRS/anti-BRS invariance and, 

s r = C FP s C FP . (17) 

On the level of the BRS and anti-BRS transformations we can have this form of Faddeev- 
Popov conjugation for arbitrary a. However, it is relatively easy to verify that the La- 
grangean, i.e. the measure of the theory, is not invariant under C PP and thus ghost and 
antighost DSEs are not identical, unless a = 1 or A = 0: With the above Faddeev-Popov 
conjugation rule for the 5-field, the sign change in the last term of (JTTj) is exactly compen- 
sated by the first term, 

iBdA + iZx(l-a)-dA(cx c) ^ i(B + Z x {\ - a)(c x c))dA - %Z X {\ - a) - dA (c x c) 

= iBdA + iZi(l-a)-dA(cx c) . (18) 

In this way, the violations of Faddeev-Popov conjugation invariance can entirely be moved 
into the term oc XB 2 , and they thus obviously disappear in the Landau gauge A = 0. On the 
other hand, in the more general ghost anti-ghost symmetric case, with a = 1 and C FP B = B, 
the theory does have the invariance under Faddeev-Popov conjugation for all A and we can 
then immediately conclude that expectation values of C FP -odd operators vanish. 
Let us now look at one of the ghost DSEs, e.g., from Eq. (fTl^) we obtain 

5 ab S xy = < |§ <* ) = Z 3 (i{dD r c) a x c b y )-^( (s r B a x ) c\ > . (19) 
oc x z 3 

For the second term on the r.h.s. we write, 
((s r B a x )cl) = (s r (BZ<*)) - (B a x (s r cD) = -{B" x B b y ) + Z x (l - |) (B« (c x c)J > , (20) 
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where we have used that expectation values of total BRS variations vanish. For the S-field 
correlations, and with its equation of motion Z\\ B = iZ%dA, one furthermore has, 

^(B a x B b y ) = 6 ab 5 xy - l^(dAldA b y ). (21) 

Inserting Eqs. (|2U|l and (|2*T|) into the ghost DSE (JT§|l we arrive at 

( dA a x dA\ ) = Z 3 ( i{dD r c)t c b y ) + iZ x (1 - |) ( (c x c )\ ) . (22) 

In the last term herein we inserted the e.o.m. for the 5-field again. This term is odd under 
Faddeev-Popov conjugation and thus vanishes in the ghost /anti-ghost symmetric case a = 1, 
as asserted above. We thus have the important form of the ghost DSE in the Faddeev-Popov 
symmetric formulation (in which there is only one such DSE), 

^{dA a x dA b y ) = Z s {i{dD r c)Z<*). (23) 

Note that we obtain the same equation for standard Faddeev-Popov theory (a = 0). The 
important difference to the standard form of the ghost DSE is given by 

Z 3 (dA a x dA b y ) - Z x X6 ab 6 xy (24) 

which vanishes in the usual Faddeev-Popov theory. For general a, however, the Slavnov- 
Taylor identities are modified also and this contribution does no longer need to vanish as we 
will see at the end of this section. Before that, we give a convenient (symmetrised) form of 
the ghost DSE valid for arbitrary a without ghost anti-ghost invariance. Note that we could 
equally have started from the ghost derivative in Eq. (JT3J) and ( c b (S/Sc x )S ) = 5 ab 5 xy . This 
would lead us to the Faddeev-Popov conjugate of Eq. (obtained from (j22j) with c — ► c, 
c — ► — c and a — > 2 — a). Adding the two, we obtain a Faddeev-Popov symmetric version in 
the place of Eq. (22J), 

A ( dA a x dA\ ) = Z 3 I (( i{dD r c)l c b y ) + ( c b yl (dD r c) a x )) - iZ x \ (1 - a) { dA a x (c x c)J > . 

(25) 

Just as we have a doubling of ghost DSEs, in absence of Faddeev-Popov conjugation invari- 
ance, we also have a doubling of Slavnov-Taylor identities. As the result of one such new 
Slavnov-Taylor identity we will derive below that 

Zi \ ( dAl (g x C )J ) = Z z % - (( i{dD r c)t <* } - ( c\ i{dD r c)l )) . (26) 

This allows us to write for the ghost DSE (|23|) and general a, finally, 

^ (9AldA b y ) = Z,{{l-^){i{dD r c)l^) + ^(c b y z(dD r cy x )) . (27) 

For a = (or 2) the l.h.s. reduces to unity and one obtains the ghost DSE of standard 
Faddeev-Popov theory. For a = l both terms on the r.h.s. are identical and add up to that 
of Eq. (USD 
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The main difference, as compared to ordinary Faddeev-Popov gauge, in an explicit rep- 
resentation of the ghost DSE will be new type of diagrams generated by the four-ghost 
interaction. The formal structure of the gluon DSE, on the other hand, remains unchanged. 

For completeness we have provided a derivation of the ghost DSE starting directly from 
the Lagrangean JI} in appendix A. For all details the interested reader is refered to this 
appendix as well as appendix B which contains the definitions of Green's functions and the 
decompositions of full into connected and one-particle irreducible Green's functions. Em- 
ploying the definitions of the bare ghost-gluon and the bare four-ghost vertex, see appendix 
B, the Dyson-Schwinger equation for the ghost propagator in coordinate space reads: 



— Z\ j d 4 zd 4 ud 4 vd 4 z x d 4 z 2 d 4 z 3 

Yf bd \y )U ,v)D%{v - Zl )Tl ha (z l7 z 3} x)D h G d (u - z 3 ). 

— Z 4 J d 4 ud 4 vT^ dfa (x,u,v,y)D f G d (v - u) 

— Z4- J d 4 zd 4 ud 4 vd 4 Uid 4 u 2 d 4 u 3 d 4 U4T^ d9 \y, z,v,u)D G e (u — U4) 

xD 9 G (v - M 2 )r^(M 3 , x, m 4 , u 2 )D J G (u 3 - z) 



^ 4 2 J ^ z ^ u ^ v ^ Ul ^ U2 ^ u ^ u ^ u ^ 



)q{V — U, 2 ;-L ^ K LOJ, L0 3 , LU 2 J-^Q 



xT k v a \u 4 ,x,u,)D^(v - u 2 )T^(u u u 3 ,u 2 )Dl d (u 3 - z). 



f28l 



Fourier transformation to momentum space yields: 
\D G {p)]- 1 = Z 3 [D§\p)]' 4 

~ Zi^f 4 f Ar?(p,?) D^ip- q)r v (q,p)D G (q) 

J d 4 q 1 q 2 T i ^ h D G {q 1 ) D G {p - q x - q 2 ) T 4gh (p, q u q 2 ) D G {q 2 
J d 4 q x q 2 r^fc D^p - q x ) D G (ql) 



9 1 9 A N 2 C 
— ^4-- 



l 4(2vr) f 

xT v (p, qx) D G (q 2 ) T^-p + q x + q 2 , ftj) D G (p -q x - 9a)- 

(29) 

The color traces have already been carried out and the reduced vertices defined in appendix 
B have been used. The four-ghost interaction generates three new diagrams in the ghost 
equation, a tadpole contribution and two two-loop diagrams. Furthermore the bare ghost- 
gluon vertex depends on the gauge parameter a, 

T^ abc (k,p,q) = gr bc (27r) 4 5 4 (k + q-p)T^(p,q) 

= [(l-£) ffM+ «pl (30) 
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FIG. 1: The coupled gluon and ghost Dyson-Schwinger equations from a BRS and Anti-BRS 
symmetric Lagrangean. Each equation contains one-loop diagrams, a tadpole contribution and a 
sunset and a squint diagram. 

Note the symmetry between the ghost momentum p p and the antighost momentum g M , when 
the gauge parameter a is set to one. 



Projection of the gluon equation 

The respective equation for the gluon propagator is formally the same as in the Faddeev- 
Popov case. Differences occur in the explicit form of the bare ghost-gluon vertex and the 
dressed vertices in general depend on the gauge parameters. The gluon DSE reads 



[D{P%1 = Z,[D^{p)\-l 

+ f d'qYf{p,q)D G {p-q)Y v {q,p)D G {q) 

1 g 4 N 2 f o 

1 g^N^ f 
- Z 4 - -— f J <? qi q 2 Tf paX D pp ,{p - qx - q 2 ) D aa ,{q 2 ) 

XlVox'O - qi - g 2 , Q2 ) D (C (p - qi ) r fvA '(p - H, qi) D x \'(qi). 

(31) 

Both equations are shown diagrammatically in fig. ^ One clearly sees the striking sim- 
ilarity between the ghost and the gluon equation once a four-ghost interaction has been 
introduced. Both equations have bare and one loop parts, a tadpole contribution, a sunset 
and a squint diagram. 
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T.I 




rsun _i_ J/sun i -f /squint . -rrsquint . -irsquint . rrsquint , -rrsquint 
TLL ' * LLL ~+~ VtTTT ~+~ TTTL ~+~ TTLL ' TLLL ~+~ LLLL 





FIG. 2: Various contributions from the respective diagrams in the transverse and longitudinal gluon 
equation and the equation for the ghost dressing function. 
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In order to sort the various contributions of the gluon equation to the inverse of the gluon 
propagator on the left hand side we project the equation on its longitudinal and transverse 
parts. It is well known that for linear covariant gauges, a = 0, the longitudinal part of the 
gluon propagator remains undressed j3|. However, away from linear covariant gauges this is 
not the case as can be seen from the corresponding Slavnov- Taylor identity derived in 
We then have three dressing functions in the general case and the propagators are given by 

DM = \D„„{p)\t + \DM\l 

= ( V _^W + AL(J)2 )M. (32) 

D G ( P ) = (33) 

The transversal and longitudinal gluon dressing functions Z(p 2 ) and L(p 2 ) can be extracted 
by contracting the gluon equation with the transversal and longitudinal projector respec- 
tively. The results are given graphically in fig. EJ where we also specify our notation for the 
different contributions being analyzed in the next section. Contributions in the transversal 
part of the gluon equation are denoted by the symbol V, contributions in the longitudinal 
part by W and the ones in the ghost equation by U. The subscripts T and L indicate the 
respective parts of the gluon propagator running around in the loops of the diagrams and 
abbreviations for the diagrams are used. For example the symbol denotes a contri- 

bution from the sunset diagram to the longitudinal gluon equation with two longitudinal 
and one transverse part of the gluon propagator running in the loop. To isolate the dressing 
functions the left hand side of the equations have already been divided by factors of 3p 2 and 
p 2 respectively. 



Generalized Slavnov Taylor Identities 



To derive the generalization of the Slavnov-Taylor identity for the longitudinal gluon 
propagator we start from the following BRS variations: 

SrpAfi) = -Z 3 {dD r c)l c\ + dAl (B-Z^icx c)) h y , (34) 
Sr(dAtc h y ) = -Z, {dD r c)lc b y - dAl {B + Z, (1 - |) (g x c)) b y . (35) 

The corresponding vacuum expectation values vanish, and taking combinations of the ex- 
pectation values of these equations we obtain: 

= (l-f)<ar(M£3)> - f (M dA °A)) 

= -Zs (1 - - ) ( (dD r cr x c b y ) + Z 3 -( {dD r c)l c\ ) + ( dA% B\ ) . (36) 

Upon insertion of the e.o.m. of the B- field, Z\\B = iZ 3 dA, this directly leads to Eq. (|27j) . 
On the other hand, the ghost DSEs from Eqs. (|13|) and (|14j) allow to eliminate the first two 
terms on the r.h.s., multiplying to them appropriate factors of | and 1 — f and inserting 
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these expression in eq. ()36|) yields 



Z 3 (dA a x dA b y ) = Z x \{s ab 6 xy (37) 
- iZ x f (1 - f ) ( IdA x c) S c b y ) + ^Z^(l- |) 2 ( ((c x c) x c): cj ) 
+ <Zx f (1 - f ) ( (M x c)% c\ > + ^ Z? ^ (1 - | ) ( ((c x c) x c): C ; ) } . 

This generalizes the Slavnov-Taylor identity for the longitudinal part of the gluon propagator 
which, contrary to the standard Faddeev-Popov gauges, does in general acquire renormal- 
ization by the interactions, c.f. Eq. (|7jl. On the r.h.s. of the Slavnov-Taylor identity, the 
terms in the 3rd line are the Faddeev-Popov conjugate of those in the 2nd. In the ghost 
anti-ghost symmetric case for a = 1 they are identical. In this case the Slavnov-Taylor 
identity simplifies, 

Zs{dA a x dA b y ) = Z x \{5 a % y - l Z 1 1 -((dAxc):4) + ^^(((cxc)xc) a x c b y )} . (38) 

Note that close to the Landau gauge the corrections to the unity of the standard Faddeev- 
Popov gauges are suppressed by one order in the gauge parameter A. 

A further Slavnov-Taylor identity is obtained by adding the expectation values of the 
BRS variations in Eqs. (|34j) and (J33j) : 

= { 8r JdAZ<*))+ (s r (dA%cl)) 

= -Z 3 ( (dD r c) a x c b y ) - Z 3 ( (dD r c) a x c\ ) - Z 1 ( dA a x (c xc){). (39) 

This leads to Eq. (|2T)j) as promised in the previous subsection. 

These Slavnov-Taylor identities indicate that the Landau gauge limit A — > is smooth. 
Based on Eqs. (}3Tj) and (}3T?j) one may anticipate that an infrared massless-like longitudinal 
part of the gluon propagator leads for sufficiently small values of the gauge parameter A to 
the same infrared enhancement of ghosts as observed in the Landau gauge. 

INFRARED ANALYSIS WITH BARE VERTICES FOR ARBITRARY GAUGE 
PARAMETERS 

In this section we will analyse the behaviour of the two-point functions at small mo- 
menta p 2 . We will employ a truncation scheme that successfully has been applied in the case 
of Landau gauge 0, IE El an d explore its applicability to general gauges. 

An interesting result of the investigations in Landau gauge is the observation, that there 
is no qualitative difference of the solutions found with bare vertices or with vertices dressed 
by the use of Slavnov-Tavlor identities. This has not only been found in truncations using 
angular approximations 0, HI for the integrals, but has been confirmed recently for a range of 
possible vertex dressings in a truncation scheme without any angular approximations . The 
reason for this somewhat surprising result has been attributed to the non-renormalization of 
the ghost-gluon vertex in Landau gauge, that is Z\ — 1. It seems as if the violation of gauge 
invariance using a bare vertex is not that severe in Landau gauge such that the resulting 
equations still provide meaningful results. In the following we will explore to what extent 
such a simple truncation idea is applicable in other gauges where Z\ ^ 1. 
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In Landau gauge the coupled set of Dyson-Schwinger equations is solved by pure power 
laws for the ghost and gluon dressing functions. Such solutions are determined analytically 
by plugging a power law ansatz in the equations and match appropriate powers on the left 
and right hand side. Once several power solutions have been found the remaining task is to 
single out the one matching the numerical solution of the renormalized equation. In Landau 
gauge it has been shown that indeed only one of the power solutions found in refs. 0, H| is the 
correct infrared limit of the renormalized solution |9j by solving the equations numerically 
for all momenta. In the following we will investigate whether there are power solutions at 
all using bare vertices for general gauge parameter a and A. 

Now we employ the power law ansatz for the dressing functions, 

G(x) = Bx 13 , Z(x) = Ax a , L(x) = Cx p , (40) 

where x = p 2 has been used. Together with the expressions for the bare vertices given in 
appendix B we plug the power laws into the ghost and the gluon equation. The formulae 
for the various integrals are given in appendix C. The straightforward but tedious algebra is 
done with the help of the algebraic manipulation program FORM |4£j • In ref . it has been 
shown that the renormalization functions Z 3 and Z 3 do not play a role in the determination of 
possible power solutions of the equations in the infrared region of momentum. Furthermore 
the tadpoles just give constant contributions to the respective propagators which vanish in 
the process of renormalization. Thus we safely omit them in the present investigation. 
For the most general gauges, a^O and A ^ 0, we obtain the following structure: 

B^x' 13 = x a+(3 (U% ress + x p+l3 (U')f' ess + x 3l3 (U') sun 

+ x a+3 yuy qumt + X p+W(uyi umt (4i) 

A- x x - a = x W(y') 9h0St + x 2a (V')$™ + x a+p (V') 9 ^ + x 2p (y') 9 ll e 

+ x 3 yvy^ T + x 2a+p (vy™ L + x a+2p {vy T u £ L + x 3p (vy™ L 

j_,r.4a (jr/\squint , 3a+p /Tr/\squint , 2a+2p f-rr/\ squint 
rx \ v )TTTT' X \ v JTTTL' x \ v ITTLL 

+x *+3 P ( V ,ywmt + x 4 P( y,yqmnt (42) 

(cxy 1 x~ p = +x 2/3 (wy host + x 2 ywy^ e + x a+p (wyy L e 

+x 3u {w')^ T + x 2a+p {wy^ L + x a+2p (wy^ L + x 3p (w')T? L 
+ x 4 ywy T q ^ + x 3r7+p (wy T q ^l + x 2u+2p {wy^y L 
+x a+3p {wy q ™y L . (43) 

Here the primed quantities are momentum independent functions of (3, a and p, c.f. fig. 121 
where the corresponding unprimed, momentum dependent quantities have been introduced. 
The pattern of the equation is such that each primed factor on the right hand side is ac- 
companied by the squared momentum x to the power of the dressing function content of 
the respective diagram. In appendix D we demonstrate how such a pattern emerges for 
example from the sunset diagram in the ghost equation, (U) sun . Note that the contributions 
(W'Yll an d (WYllll are zero an d therefore missing in the longitudinal gluon equation (f4*3l 
as momentum conservation cannot hold with three longitudinal gluons in the three gluon 
vertex. 

For the following argument we focus on one particular contribution on each right hand 
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side of the equations: 



B^x-P = x w {U') sun + ... (44) 

A~ x x~ a = x 4ct (K')tttt + • • • ( 45 ) 
(CXy'x-P = x i »{W')T£ L + ... (46) 

The coefficients (U') sun , (V')tttt an d (W'Ylll are nonzero and explicitly given in ap- 
pendix D. First it is now easy to see from equations (|44jh ()45|) and ()46|) that neither (3 nor 
a nor p can be negative. If one of these powers would be negative the limit x — > would 
lead to a vanishing left hand side of the respective equation whereas the right hand side 
is singular in this limit. This is a contradiction as the power on the left hand side of the 
equation should match the leading power on the right hand side. Second if one of (3, a or p 
would be positive, then the diverging left hand side of the respective equation would require 
a diverging counterpart on the right hand side. However, all powers on the right hand side 
are positive as we already concluded that (3, a or p are not negative and there are no minus 
sines in any powers on the right hand sides, c.f. Eqs. (JUJ), (J4*2*|) and Therefore for 

positive powers all terms on the right hand side vanish in the limit x — > which leads again 
to a contradiction. The last possibility is then (3 = a = p = 0, but then one gets pertur- 
bative logarithms on the right hand side of the equation which do not match the constant 
on the left hand side. Thus in the all-bare-vertex truncation there is no power solution for 
general values of the gauge parameters A ^ and a^O. Based on the considerations on the 
Slavnov-Taylor identities given in the previous section we therefore arrive at the conclusion 
that this truncation is insufficient to determine the infrared behaviour of the propagators 
even qualitatively. 

There are two limits for the gauge parameters a and A in which the situation changes. The 
first one is a = 0, that is ordinary linear covariant gauges. Due to the corresponding Slavnov- 
Taylor identity the longitudinal part of the gluon propagator remains undressed, L(p 2 ) = 1 
However, replacing dressed vertices by bare ones in the infrared, this identity might be 
violated (which does not happen in perturbation theory, of course). We therefore employ the 
general expression L(p 2 ) = C(p 2 ) p for the longitudinal gluon dressing function and explore 
whether the limit p — > can be taken with bare vertices. In the ghost equation the squint as 
well as the sunset diagram disappear and we are left with the one-loop contributions U^ ress 
and U[ ress . The explicit expression for the ghost equation is given by (cf. appendix D) 

B^x'P = u^ ress + Ut ress (47) 
= x(0 ^ g 2 N c Z x AB -3 T(2 + g) T(l + a) T(2 - (3 - a) 



16tt 2 2(J3 + <t){-1+/3 + a) Y{1 - (3) T(2 -a)T(3 + (3 + a) 
( ^ p) g 2 N c Z x \BC p + 1/2 T(2 + P)T(l + p)T(-P-p) 

16tt 2 T(-P)T(2-p)T(3 + p + P) ' 1 ' 

For p —>■ we run into the same contradiction as explained above for general values of the 
gauge parameters a and A. However, admitting the generation of a (spurious) longitudinal 
gluon dressing this contradiction can be resolved in the following way: Equation (J45j) for the 
transversal gluon dressing function does not change in structure, therefore o > 0. Further- 
more we have p > from eq. (J46j) . Then we have —2/5 = a and/or —2/5 = p in the ghost 
equation ()48j) and P < 0, i.e. a diverging ghost dressing function in the infrared. From this it 
follows immediately, that the ghost loop is the dominant contribution in both, the equations 
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for the transversal and longitudinal gluon dressing function. From these two equations we 
therefore infer 

-P = a/2 = p/2 =: k, (49) 

which is consistent with the ghost equation. We thus find an infrared vanishing gluon dressing 
and a singular ghost dressing function for all values of the gauge parameter A. This result 
is identical to the one in Landau gauge @, E| . However, a word of caution is in order. In 
Landau gauge there are indications |4| |41( that the general result ()49|) does not change when 
the vertices are dressed. This has been confirmed recently for a range of possible vertex 
dressings (§[. It is an as yet open question whether this is true for A 7^ in the same way. 

Having addressed the case of linear covariant gauges with a = we now turn to the other 
interesting limit, that is A = 0, while a / 0. It is easy to see, that the a-dependence of 
the Lagrangean can be eliminated in this case by partial integration using the constraint 
dA = 0. However, on the level of the DSEs with bare vertices there remain spurious in- 
dependent terms on the right hand side of the gluon equation. In the next section we will 
investigate the dependence of the Landau gauge solution on these spurious a-terms. 



SOLUTIONS IN LANDAU GAUGE 

To assess the influence of the spurious a-terms in Landau gauge we use the truncation 
scheme developed in |9[. There the two loop diagrams in the gluon equation have been ne- 
glected as they are subleading in the perturbative regime and ghost loop dominance has been 
assumed in the infrared. In order to obtain the correct one loop behaviour of the ghost and 
gluon dressing functions the gluon loop has been modified by replacing the renormalization 
constant Z\ by a momentum dependent function Z\: 

Z X (L, s) -> Z x (x, y, z- s, L) = Z ^ )(1+3S) Z{z)( i +3S) ( 50 ) 

Here L = A 2 denotes a cutoff and s = fi 2 a renormalization scale in units of squared momenta. 
The momentum x = p 2 is the one flowing into the loop, y = q 2 is the loop momentum over 
which is integrated and z '■— k 2 = (p — q) 2 . Furthermore the anomalous dimension 5 of the 
ghost dressing function has been used. The gluon equation is contracted with the general 
tensor 

P#(P) = ^-C^- (51) 

As a completely transversal gluon equation would be independent of the parameter ( the 
use of the general projector provides an opportunity to test for violations of transversality 
due to the truncation. For ( 7^ 4 one has to take care of spurious quadratic divergencies that 
have to be subtracted in the kernel of the gluon equation. 

The coupled set of equations for the ghost and gluon dressing functions then read as 
follows 

1 - W*. / A^^GMZM , (52) 



G(x 



(27T)* 


xy 


d 4 q 


M(x,y,z) { 


(2tt) 4 


xy 


f d 4 q 


Q(x,y,z) 



2 N C f d A q Q{x,y,z) 

+g T / 7^4 Z{y)Z{z)Z 1 {y,z) . 53 

3 J (2tt) 4 xy 
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FIG. 3: Here the graphical solution to equation (|58|) is shown. The thick line represents the left 
hand side of equation (|55|) , whereas the other curves depict the right hand side for different values 
of the parameters The left figure shows results for a = and a = 2, whereas in the figure on 
the right a = 1. The ellipse marks the bulk of solutions between k = 0.5 and k = 0.6 for £ = 1, 
whereas the circles in the left figure show the movement of the solution for the Brown-Pennington 
case C = 4 from k = 1 to k = 1.3. 

The kernels ordered with respect to powers of z := p 2 = (k — q) 2 have the form: 



K(*>V>*) = -A —) + z{~) "4 (54) 



M(x,y,z) = - x H H 55 

v ' y ' ; z V 4 2 4 x / 2 2x Ax V ; 

i /h 3 2 19- c 5-c 2 (y 3 

Q{x,y,z) = — [-— + x — xy + — — y + - — 

z l \ 8 y 8 4 8s 



l fx 2 15 + C 17-C M ,y 2 

+- — x — y + C— 

z \ y 4 4 x 

19-C s , 17-C ] %y 

8 y 4 4 x 

(~ + ^) + z 2 ^~ + 7(4 - 0- (56) 
\x Ay J 8xy 4 V 7 v 7 



First we accomplish the infrared analysis. With equation (}4T?j) we employ the ansatz 

Z(x) = Ax 2K G{x) = Bx~ K (57) 
in the equations (|52|) and (|53p . After integration we match coefficients of equal powers on 
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both side of the equations and obtain 

1 (2 + k)(1 + k) 
18 (3 - 2k) 

(4k -2) (-1 + k) 



(C - 1) [4 k 2 (a 2 - 2 a + 1) + 8 k a (2 - a) + 3 a (a - 2)] + k (10 - 7 C) - 6 + 3 C ' 

(58) 

The values of k for different projectors V^ 1 can be read off fig. El The curve given by the fully 
drawn line represents the term on the left hand side of equation (J58)) . whereas the other lines 
depict the right hand side for several values of the parameter (. Only the two ( = 1 solutions 
are manifestly independent of a, as pointed out in (if. The spurious a-dependence of the 
( = 4 values reported therein, here implies that general ( solutions must necessarily show 
such an a-dependence also, whenever ( ^ 1. However, the bulk of solutions between k = 0.5 
and k = 0.6 remains nearly unchanged when a is varied, whereas most of the solutions for 
k > 1 disappear. For the Brown-Pennington projector ( = 4 no solution can be found for 
the symmetric case, a = 1, in complete agreement with the findings of ref. Indeed it has 
been shown that only the smaller solutions are those that connect to numerical results 
for finite momenta. 

We now explore the impact of the spurious a term on the behaviour of the solutions for 
all momenta x. We have solved eqs. (p)2*|) and (p)3*|) numerically using the same technique 
as described in The results can be seen in fig. 0J As the dependence of the kernel of 
the ghost loop on a vanishes in the case of the transverse projector, ( = 1, this solution 
is the same as the one already calculated in For the other cases the power k changes 
from 0.5953 for £ = 1 to 0.5020 for £ = 3.9 in accordance with the infrared analysis. The 
ultraviolet properties of the solutions are slightly disturbed compared to the cases a = 
and a = 2. An analysis of the ultraviolet behaviour done similarly to the one in ref. [j| 
reveals that the a-term in the ghost loop induces a spurious dependence of the anomalous 
dimensions on the parameter £: 
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-26 - (C - l)of(2 - a) 
44+ (C- l)a(2-a) 
-9 

(59) 



44+ (C - l)a(2-a)' 



For general a only the transverse projector removes the spurious term in the ghost equation 
and leads to the correct one loop scaling of the equations, that is S = —9/44 for the ghost 
and 7 = —13/22 for the gluon dressing function for an arbitrary number of colors and zero 
flavours. 



CONCLUSION 

We have studied the infrared behaviour of the ghost and gluon propagators in general 
covariant gauges. These gauges allow to interpolate via a second gauge parameter between 
the linear-covariant ones of standard Faddeev-Popov theory and the ghost-antighost sym- 
metric gauges. We derived the corresponding generalised Dyson-Schwinger equations for 
the propagators which include the ones of linear-covariant gauges as the limit where the 
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G(x) 
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x=p 2 [GeV 2 ] 

FIG. 4: Shown are the gluon dressing function, the ghost dressing function and the running 
coupling in the truncation scheme 0| for the gauge parameters a = 1 and A = and different 
projectors V^'- 

second gauge parameter vanishes. Note that ghost-antighost symmetric gauges are particu- 
larly interesting as they allow an interpretation of the antighost field being the antiparticle 
of the ghost which includes also the possibility of ghost-antighost condensate. Due to the 
emergence of a four-ghost interaction term in the Lagrangean for general values of gauge 
parameters the Dyson-Schwinger equation of the ghost propagator displays a rich structure 
very similar to the one of the gluon equation. On the other hand, in the gluon equation we 
obtain the same structure as in linear covariant gauges apart from the fact that the gluon 
propagator acquires a nontrivial longitudinal part which appears in turn in all diagrams. 
The gluon and ghost equations depend therefore on three dressing functions, one for the 
ghost, one for the transverse part of the gluon propagator and one for the longitudinal one, 
which are constrained, however, by Slavnov- Taylor identities in an intricate way. 

We then employed a truncation scheme for the Dyson-Schwinger equations that uses bare 
vertices in place of the dressed ones. The success of this particular truncation scheme in 
Landau gauge has been attributed to the non-renormalization of the ghost-gluon vertex, 
that is Z\ = 1. We addressed the infrared behaviour of the ghost and gluon propagators 
for general gauges by employing power law ansdtze for the respective dressing functions. 
We then have been able to evaluate the infrared behavior of the gluon and ghost equations 
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analytically. 

For all linear covariant gauges we find a similar result as compared to the one in Landau 
gauge: An infrared suppressed gluon propagator and an infrared enhanced ghost. Whereas in 
Landau gauge there are indications that this generic result is not changed when the vertices 
are dressed |8J], it remains an open question whether this is the case in linear covariant gauges 
in general. Away from linear covariant gauges, that is in the general case a ^ and A 7^ 0, 
we do not find power solutions for the dressing functions. However, we expect this to change 
with appropriate vertex dressings. Nevertheless, it remains to be emphasized that therefore 
also the occurrence of a ghost and/or gluon mass is excluded in this specific truncation 
scheme within this class of gauges. A Dyson-Schwinger equation based investigation of 
the related question of a ghost-antighost vacuum condensate, or more generally, of an 'on- 
shell'-BRS-invariant dimension two condensate, needs to take into account the generalized 
Slavnov-Taylor identities (|37)l and ()39|). The question arises whether an infrared massless- 
like longitudinal part of the gluon propagator leads for all values of the gauge parameters 
to the same infrared enhancement of ghosts as observed in the Landau gauge. Work in this 
direction is in progress. 

A special case among all gauges considered here is Landau gauge. In the limit A = 
the general Lagrangean ((TJ) becomes independent of the second gauge parameter a, thus 
Landau gauge is also a special case of ghost-antighost symmetric gauges. Although the 
Lagrangean of the theory is independent of the gauge parameter a, our simple truncation 
scheme breaks this invariance and spurious a-dependent terms arise in the ghost loop of the 
gluon Dyson-Schwinger equation. Examining the case a = 1 we showed that the influence 
of these spurious terms is very small. We determined solutions for the ghost and gluon 
dressing functions both analytically in the infrared and numerically for finite momenta and 
found solutions identical to the ones of ref. 0| provided the gluon equation is projected onto 
its physical, transversal components. We thus recovered the results of Landau gauge from a 
different direction in the two dimensional space of gauge parameters. 
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APPENDIX A: DERIVATION OF THE DYSON SCHWINGER EQUATION 
FOR THE GHOST PROPAGATOR 



We start by transforming the Lagrangean ((TJ into a more suitable form by partial integra- 
tion, assuming the usual boundary conditions of vanishing fields at infinity. In order to keep 
notation on a readable level we will suppress renormalization constants in this appendix: 
The derivation of the Dyson-Schwinger equation for the ghost propagator remains formally 
unchanged by the rescaling (jSJ) and thus the appropriate renormalization constants can be 
regained straightforwardly. We obtain 

£ = l -Al (-d 2 6» u + (l - £) d,d)j A a u - gf abc (d,Al) A\A C V 

+^f abe f cde A»A b v A*Ai + c a d 2 c a + |(l - I) -g 2 r e f bde c a c b c c c d 
+i (l - |) gf abc c a d, (A c /) + i^gr hc c a A^d»c h . (60) 
The partition function of the theory is given by 

Z[J, a,a} = J V[Acc] exp j- J d A zC + J d A z (A a J a + ac + ca) J (61) 

with the sources J, a and a of the gluon, antighost and ghost fields, respectively. The action 
is given by S[J,c,c] = J d^zL. The generating functional of connected Green's functions, 
W[J, a, a], is defined as the logarithm of the partition function. The functional Legendre 
transform of W is the effective action 

T[A, c, c] = -W[J, cr,a} + J d 4 z {A a J a + ac + ca) , (62) 

which is the generating functional of one-particle irreducible vertex functions. The fields and 
sources can be written as functional derivatives of the respective generating functionals in 
the following way 

5W _ 5VV_ 5W 

5a 5a 5J^ M ' 

Tc =a ' Tc =a > Ja, = j - (63) 

The sign conventions have been chosen such that derivatives with respect to c and a are left 
derivatives whereas the ones with respect to c and a are right derivatives, 

5 5 

left derivative —/ \ := right derivative. (64) 



5 (cr, c) 5 (a, c) 

Given that the functional integral is well-defined, the Dyson-Schwinger equation for the 
ghost propagator is derived from the observation that the integral of a total derivative 
vanishes provided the measure is invariant under field translations. We take the derivative 
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with respect to the antighost field and obtain 



ac + ca , 



= J V[Acc]-^exp j- J d 4 zC + j d 4 z (A a J a + 

= j V[Acc] (- — [ ^ C,C] + a\ exp | - J d 4 z £ + J d 4 z (A a J a + ac + ca) J 

[sj ' s *' 5 °±+a)z[J,a,<T]. (65) 



5c 

Now we use the relations (jBSj) and apply a further functional derivative with respect to the 
source a b (y). We arrive at 

with explicit colour indices and space-time arguments. Setting the sources equal to zero we 
obtain the ghost Dyson- Schwinger equation 

!)S c b {y)\ = 5{z-y)5 cb . (67) 



5c c [z i 

The derivative is easily calculated 

= <9V(z) + f (l - f ) \g 2 r de i he c\z)c^z)c\z) 

+t i 1 ~ f ) 9fCded » ( A W^ Z )) + ^9f C ^(z)d,c d (z). (68) 

Whereas in the covariant formalism full and connected three-point functions are the same, 
the four-point correlations have to be decomposed into disconnected and connected parts. 
For the four-ghost correlation function this results in 

(c b (y)c d ( z y(zy(z)) = (c*{ v y{ z )) (<?*(*)</(*)> - (c b ( y y( z )) (c d ( z y(z)) 

+ (c b (y)c d (z)c f (zy(z)) conn .. (69) 
Keeping in mind the Grassmann nature of the ghost and antighost fields we then obtain 

- 6(z - y)5 cb = d 2 (c b (y)c c (z)) + | (l - | ) ^g 2 f cde f^ {(<*(y)<*{z) C f(z)d'{z)) 
+ ((c b (y)c°(z)) (c d (z)c^z)) - (c b (y)c^(z)) (c d (z)c° (z)))} 

)d,c d {z)) , 

(70) 



+ (l - f ) 9f cde (c\y)d, {A^(z)c d (z))) + ^9f cde ^(y)A^z)d,c d (z)) 



where all correlations are connected Green's functions. We now use the relation 



5a a (x) J 5c d (z) 5a a (x) J 5c d (z)5c b (y) 5a a (x)5a d (z) 

J d 4 z [D db {z - y)] ^ Dg(x - z) 

(71) 
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and multiply eq. ([70]) with — [Dq(x — z)\ 1 = [(c c (z)c a (x))} \ We arrive at 
[Dgix-y)}- 1 = d 2 5(x-y)5 ab 

+ (c b (y)c°(z)) (c d (z)c f (z)) - (c\y)c f (z)) (c d (z)c°(z))} 

- ^( 1 -f)^ e / ^[^-^)]- 1 (c 6 (y)^(^(z)c d (z))) 

- i\gf cde I d*z[D a G ?(x - z^ic^A^d^iz)). 

(72) 

Before we decompose the connected Green's functions into one particle irreducible ones we 
have to take care of the space-time derivatives. Noting that 

5 2 W f „ 8 2 W 



»5Jc(z)o d (z) J ^ vv " 5Jc(u)a d (u) 



- I d*iuv]d; (8(u - z)) 8{u - v) 5Jc ^ d(u) (73) 



with the abbreviation d 4 ud 4 v =: d 4 [m>], and 



8 8W f 4 5 au 8W 

: d l-rrr = d [MS(u-z)8(u-v)——rd^- 



5Jc{z) **a d (z) J ' y 8J*(v) "a d (u) 



= -j d 4 [uv] d; (8(u - z)8{u - v)) 5J !^ ad{u) (74) 

we can replace the derivative terms by the bare ghost-gluon vertex defined in appendix B. 
The tadpole term can be treated in the following way: 

*z[Dg(x - z)]- 1 r de f^{(c b ( y y(z)) (c d (z)cf(z)) - (J(y)cf(z)) (c d (z)c°(z))} 
= 2 J d 4 z [D a G c (x - z)}- l f cde f f9e {(c b (y)c 9 (z)) (c d (z)c f '(*))} 
= 2 J d 4 [zuv] [D a G c (x - z)}-^^ - u) 8{u - v) f cde f^ e D 9 G (z - y) D f G (v - u) 

= 2 J d 4 [uv] 5(x - y) 5(z - u) 8{u - v) f bde ff ae D f G d {v - u). (75) 

Plugging the expressions for the ghost-gluon loop and the one for the tadpole into eq. (J72jl 
and using the expression for the bare four-ghost vertex given in appendix B we obtain 



[Dg^-vT 1 = d 2 5(x - y)5 ab - J d 4 [uv]T 



Agh 0, u, v, y) D f G (v - u) 



+ f(!-f)^Wx 



d 4 [zuv] 5{z - u) 5{u - v) [D a G c {x - z)}' 1 (c b {y)c d (z)c f (u)c 9 (v)) 

- I rf 4 ^Hri 0)cde (^w^)[^ c (x-^)]-Mc 6 (y)^(^)c d H). 

(76) 
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To decompose the connected Green's functions into one-particle irreducible ones we use 
the relations 



(A;(v)c b (y)c d (u)} = J d\z x z 2 z 3 \D%{v-z x )D h £{y-z 2 )T{ h \z x ,z 3 ,z 2 )D\ 

(77) 

(c b (y)c d (z)c f (u)c a (v)) = J d 4 [ Ul u 2 u 3 u 4 u 5 u 6 ] [Df v { Ul - u 4 ) D§(u - u 5 ) 

x r* w (u 4 , u e , u 5 ) D hb {u & - y) D%(v - u 2 ) 
xYf{u u u 3 ,u 2 )D id {u 3 -z)^ 

- j d 4 [uiu 2 u 3 u i us\ ^D^iu-u^D^iu^-y) 

x D 9 q(v - u 2 ) Tl h g e h \u 3 ,u 5 ,u 4 ,u 2 ) D ]d (u 3 - z)} , (78) 

which have been derived in appendix B. 

Substituting these expressions into eq. (|7S|) we arrive at the final expression for the ghost 
Dyson- Schwinger equation in coordinate space: 

[D^x-y)]' 1 = [DP^ix-y)]- 1 

- J d 4 [uv]Tf^ dfa (x,u,v,y)D G f(v - u) 

- - I d^zuvuxu^UiUs] rf^ 9f (y, z, v, u) Df v {ui - u 4 ) D§(u - u 5 ) 



2 J [ 1 ^ o 4 5J A.fjh \y '1 A '1 v '1 u 7 " L> >iA 

x Y k v a \u A , x, u 5 ) D%{v - u 2 ) Vf (u h u 3 , u 2 ) & d (u 3 - z) 

- -j d 4, [zuvu 1 u 2 U3Ui]Tf^ f (y,z,v,u)D^(u-U4) 

x D 9 q(v - u 2 ) Yf^{u 3 , x, u 4 , u 2 ) D J c(u 3 - z) 

- J d%uvz 1 z 2 z 3 ]Yf bd \y^v)D%{v-z l )Yl h \z l ,z 3 ,x)D h G d {u-z 3 ) 

(79) 

where an additional minus signs arises from the interchange of the colour indices / and g in 
the bare four-ghost vertices and from the interchange of j and i in the ghost-gluon vertex. 
After performing a Fourier transformation we obtain the respective expression in momen- 



24 



^ a d c 

§ ^ h ''' 

§ V 

( l /\ k 2 

P ^ / kl \ 

b c ' 

a b 

FIG. 5: Momentum routing for the tree level ghost-gluon and four-ghost vertices. 

turn space 
[Dcip)}- 1 = [D^ipT 1 

+ ^)^/Ar£i£ G (g) 

+ (^) \j£y s J d*[ qi q 2 ]T^ h D, v (p- qi )D G {ql) 

x I\,(p, q{) D G (q 2 ) r M (-p + q x + q 2 , q%) D G (p - q x - q 2 ) 

- (-^c) 2(2^)8 _/ rf4 ^ 1 ^] r i i jD G(gi)£>G(p-gi-g2)r 4s / l (p,gi,g2) J DG(g2) 

+ ("^ (1^)4 / d*qTf(p,q)D, lv (p-q)T u {q,P)DG{q) 

(80) 

where the colour traces have been carried out and the reduced vertices defined in appendix 
B have been used. 

APPENDIX B: DEFINITIONS AND DECOMPOSITIONS 
Ghost and gluon propagators: 

The full ghost and gluon propagators in coordinate space are defined to be 

SW 

SW 

{A > )AM) = jjmw) =Klx - y) - (82) 

(83) 

The bare propagators in coordinate space can be easily derived from the quadratic part of 
the action, 

S quad = J d A x> (-d 2 5,„ + (i - x) 9 A) K + c a 9 2 c a | (84) 
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and are given by 
[D^ a \x - y)} 



5 2 S, 



quad 



8 ab d 2 5{x - y) 



5c a (x)c b (y) 



(85) 



-i 



5 2 S, 



quad 



5A-(x)A b (y) 



S ab [ -<9 2 <V + ( 1 - j- ) dpd v ) S(x - y) , (86) 



with the gauge parameter A. After Fourier transformation one obtains the corresponding 
expressions in momentum space: 



Df a \p) 



-i 



abJ2 



-5 ao p 



(87) 



Ghost-gluon vertex: 

From the ghost gluon part of the action 



s ghgl = J dV {-i (i - 1) <?r bc (^c) a c / + ^/"^^a"^} 



(89) 



the tree level ghost gluon vertex T abc is easily derived: 



r^ abc (x,y,z) 



6 3 S, 



ghgl 



5A a (x)5c b (y)Sc c (z) 



9f abc [i (l - f ) (^ 4 (z - y)) 5 4 (z -x) + ifa (8\z - y)5\z - x)) 



(90) 

Using the momentum conventions of Fig. the Fourier transformed bare ghost-gluon vertex 
reads 



r^ abc (k,p,q) = I d\xyz\ Tf c (x,y,z)e^ x+ ^-^ 
= gf abc (2vr)V (k + q-p) 



1 - ^ J % + 2^. 



(91) 



where the abbreviation d 4 x d A y d A z =: d 4 [xyz] has been introduced. Note the symmetry of 
the vertex in the ghost momenta p^ and if a — 1. For convenience we define a reduced 
vertex function (p, q) by 



T^ abc (k,p,q) = gf abc (2n)^(k + q-p)T^(p,q) 
rf (P,q) = 



a\ a 
1 - 2 ) % + 2^. 

The full one particle irreducible ghost gluon vertex in coordinate space is given by 

rf c (x,y,z) <)l 



5J^(x)5c b (y)5c c (z)' 



(92) 



(93) 
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Four-ghost vertex: 

The four-ghost vertex is derived from the four ghost part of the action 

" N X g 2 f ace f bde c a c b c c c d \ (94) 



Sigh z 

which leads to 

rL [x,y,z,w) 



d A x' 



a 



1 - 



2 V 2/2' 



5c a (x)5c b (y)Sc c (z)Sc d (w) 
a 



- 1) W-f-^S^&Hx - y)s'(y - z)s\z - w). 



2V 2J* - - < 95 > 

Again using the momentum conventions of figure (0) one obtains for the Fourier-transformed 
bare four-ghost vertex 

5 4 S, 



r^ abcd (h, fa, fa, fa) 



Hgh 



5c a (x)Sc b (y)Sc c (z)Sc d (w) 
= f (l-f) V/^/^W*! + k 2 - fa - fa). 
We define a reduced vertex function vf^ by 

rf g )abcd (k u fa,fa,fa) = g*f^f^(2K) 4 5\k l + fa-fa-fa)Tf g 



(96) 



r (o) 



The full four-ghost vertex in coordinate space is formally given by 

5Y 



^abcd. 



x,y,z) 



5c a (x)5c b (y)5c c (z)5c d (y) ' 

Decomposition of connected ghost-gluon Green's function: 

With the help of the matrix relation 

for 1 ^ _ -ifbc -i 



(97) 
(98) 



(99) 



d A z- 



8 2 T 



5 2 W 



and the identity 

S(v — x)5' 

J " ~ 5c d (z) 5a a (x) J " " 5c d (z)5c b {y) 5a a (x)5a d (z) 

we decompose the connected ghost-gluon correlation function, (y4^(x)c b (|/)c c (z)), in the fol- 
lowing way: 



(100) 



(A a (x)c b (y)c c (z)) 



8 



5 3 W 



5J-(x) 



5J" L (x)5a b (y)5a c (z) 



5 2 T 



5c b (y)5c c (z) 



, SAt( Ul ) 5 



5J a Ax) 5A d (u x ) 



5 2 T 



d 4 [u 1 u 2 u 3 ] 



8 2 W 



S^(y)8c c (z) 
S 2 W 



5 3 T 



5 2 W 



5Jf l {x)5J d {ui) 5a b (y)5cr e (u 2 ) dA^u^S^iu^dcf (u 3 ) 5af (u 3 )5a c (z) 



J d^u^Df^x - Ul )D e G b (u 2 - y)r d ^( Ul ,u 2 ,u 3 )D c G f (u 3 - z). 



(101) 



27 



Here we used the abbreviation d A [u\U2U 3 } := d A u\ d A U2 d A u 3 and the definitions of the gluon 
propagator D^, the ghost propagator Dq and the ghost-gluon vertex T u given in previous 
subsections. 

Decomposition of connected four ghost Green's function: 

Furthermore we need the decomposition of the four-ghost correlation function into one- 
particle irreducible parts. We start at a stage where the sources are still present and set 
them to zero at the end of the derivation. We first give the decomposition of the connected 
ghost-antighost-ghost three-point function 



{c b {y)c c {z)c d {w)) 
8 



8 3 W 



8a b {y) 
d A U\ 



8a b (y)8a c (z)8a d (w) 



8 2 T 



_8a c (z)8a d (w) 

SA e v ( Ul ) 8 
8a b (y) 5A*{u x ) 
8 2 W 



8 2 Y 



-i -i 



d [uiu 2 u 3 ] 



8a c (z)8a d (w) 
5 2 W 



8 3 T 



5 2 W 



8a b {y)8J^(u 1 ) 8a c (z)8af (u 2 ) 8Al{ux)8cf{u 2 )8c.9{u 3 ) 8<j9{u 3 )Sa d (w) ' 

(102) 



Then we decompose the connected four-ghost Green's function: 

S A W 



{c a {x)c?{y)c c {z)c?{w)) 



8a a (x)8a b (y)8a c (z)8a d (w) 



6 



Sa a (x) 



d A [uiu 2 u 3 ] 



5 2 W 



5 2 W 



5a b (y)5J £ Ju 1 



8 3 T 



5a c (z)5af (1*2) 
5 2 W 



5A e (ui)5cf(u 2 )8&(u 3 ) 8as(u 3 )8a d (w) 



(103) 



Carrying out the remaining derivative gives four terms. The two terms where the derivative 
acts on the second and on the last propagator vanish, because the term g^F^^jwn^) vanishes 
when the sources are set to zero. The contribution where the derivative acts on the first 
propagator can be treated using eq. ()101|) . In the expression with the derivative acting on 
the vertex we use 



5 2 W 



5 A T 



5 A T 



5a b (y)5J e Aui) Sa a (x)8A e (ui)Scf{u 2 )8c9(u 3 ) 5a a (x)5a b (y)5cf (u 2 )5&(u 3 ) 



d A Ui 



5 2 W 



5 A T 



d [M4M5] 



Sa a (x)a e (u 4 ) 5c e (u A )8a b (y)8ct '(u 2 )S&(u 3 ) 
8 2 W 8 2 W 



8 A T 



d [M4M5] 



8a a (x)a e (u 4 ) 8a b (y)8a h (u 5 ) 8c e (u 4 )8c h (u 5 )8c^ \u 2 )8&(u 3 ) 
8 2 W 8 2 W 8 A T 



8a a (x)(T e {u 4 ) 8a h (u 5 )8a b (y) 8c?{u 3 )8c h {u 5 )8(f{u 4 )8cf{u2)' 

(104) 
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Collecting all this together we arrive at 

f ( ^ w 8 2 W 

(c a (x)c (y)c c (z)c (w)) = / d [uiU 2 U 3 U 4 UhU&] < f r ( . 

5 3 T 5 2 W 5 2 W 

X 5Al(u 4 )Sc9{u 5 )Sc h {u 6 ) Sa h (u 6 )5a b (y) 5tT c (^)5a i (u 2 ) 

^(wi^c^Ma)^'^) 5aJ(M 3 )5a d («;) J 

- /^[«i^3«4« 5 ]{ 5<rfl( ^ M ^J^y 

5a c (z) 5a-/ (« 2 ) 5c» (tt 3 ) fc^ (u 5 ) Sc e (u 4 )6cf (u 2 ) 
5 2 W ^[ 

x ~ s 7 \ r- T/ \ ^ • (105) 
5o-9(M 3 )5a d (u7) J v 7 

Interchanging some Grassmann fields in the correlations and using the definitions for the 
propagators and vertices given in the previous subsections we arrive at 



(c\y)c d (w)c a (x)c c (z)) = j d 4 [ Ul u 2 u 3 u 4 u 5 u 6 } {D e J v ( Ul - u A ) D^(x - u 5 ) 

xTl h %u 4 , u e , u 5 )D h G b (u 6 - y) D%{z - u 2 ) 
xrf( Ul ,u 3 ,u 2 )D> G d (u 3 -w)} 

+ J d 4 [u 1 u 2 u 3 u i u 5 ] {D a (?{x - u 4 ) D G b (u 5 - y) 

xD c J(z - u 2 ) Tl 9 g e h f (u 5 ,u 3 ,u 4 ,u 2 ) D 9 G (u 3 - w)j , (106) 

which is the decomposition of the four-ghost correlation used in appendix A. 



APPENDIX C: TENSOR INTEGRALS 



The explicit expression for the scalar bubble integral /, defined in eq. f|107|) . can be easily 
evaluated in Euclidean space-time using the Feynman-parameterisation. With the squared 
momenta x = p 2 , y = q 2 and z — (p — q) 2 the result is given by 

I(a,b,p) := [ d 4 q-^- b (107) 

J y a z b 



2 2-a- b r(2 - a) r(2 - b) r(o + b - 2) 



= * - x T77 hrkTTTTi -■ (108) 

F(a) r(6) r(4 - a - b) v ; 

The corresponding tensor integrals can be reduced to scalar integrals by extracting combi- 
nations of momenta p^ and the symmetric tensor according to the symmetry properties 
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of the integrand: 
K^ u (a, b,p) 



jA 1p 
a q 



y a z b 



y a z b 



Ji(o,6,p) pp, 

Kx(a,b,p) p^p u + K 2 (a,b,p) xS^, 



(109) 
(110) 



J pvp 



M, 



pupa 



a, b,p) 



a,b,p) 



d 4 q^QuQp =^6^) p,jPuPp 

y a z° 

+ L 2 (a,b,p) x {pp5 up + Pu& P p+ Ppbpu) , (111) 

A ^^^ = Mi(a,&,p) PpPuPpPa 

+ M 2 (a, x (5p U p p p a + 5 fip p u p a + Sp a p P Pp+ 

^up PpPa + &uu PpPp + Sp a PpP v ) 

+ M 3 (a, b,p) x 2 (5p U Spa + 5p P b va + 6^ 8 pv ) . (112) 

The scalar integrals in these expressions are calculated by contracting them with appropriate 
tensors, writing all scalar products in terms of squared momenta x, y and z and applying 
eq. ()108j) . One arrives at 



Ji 
Kx 
K 2 
Ll 
L 2 
Mi 
M 2 
M 3 



7T 



7T 



7T 



7T 



7T 



7T 



7T 



7T~ 



2 r(3- 


a 


r(2-&)r< 


a 


+ b- 


2) 


r 


[a] 


r(6)r(5- 


a 


~b) 




2 r(4- 


a] 


r(2-6)T( 


a 


+ b- 


2) 


r 


a 


r(6) r(6 - 


a 


~b) 




2 r(3- 


a) 


r(3-6)r< 


a 


+ b- 


2) 


r 


a 


r(6) r(6 - 


a 


~b) 




2 r(5- 


a 


r(2-6)T( 


a 


+ &- 


2) 


r 


[a 


r(6)r(7- 


a 


~b) 




2 r(4- 


a] 


r(3-6)r< 


a 


+ b- 


2) 


r 


a 


r(6)r(7- 


a 


-b) 




2 r(6- 


a) 


r(2-6)r< 


a 


+ b- 


2) 


r 


[a] 


r(6) r(8 - 


a 


~b) 




2 r(5- 


a) 


r(3-6)r< 


a 


+ b- 


2) 


r 


a 


r(6) r(8 - 


a 


~b) 




2 r(4- 


a 


r(4-6)T( 


a 


+ &- 


2) 


r 


a) 


r(6)r(8- 


a 


~b) 





X 



2-o- 



X 



2-a-b 



2-a-b 



„2-a-b 



„2-a-& 



2-a- 



X 



X 



2-a-b 



■ X 



2-o- 



(113) 
(114) 
(115) 
(116) 
(117) 
(118) 
(119) 
(120) 



APPENDIX D: EXPRESSIONS FOR SOME DIAGRAMS IN BARE VERTEX 
APPROXIMATION 

In this appendix we give explicitly the expressions for some diagrams needed for our 
investigation in the main body of the paper. All algebraic manipulations have been done 
using the program FORM ji^ . Our ansatze for the small momentum behaviour of the ghost 
dressing function G, the transversal gluon dressing function Z and the longitudinal gluon 
dressing function L are the power laws 

G(x) = Bx 13 , Z(x) = Ax\ L(x) = Cx p , (121) 
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FIG. 6: Momentum routing for the sunset and for the dressing diagram in the ghost equation. 



where we have used the abbreviation x = p 2 . 

We first evaluate the sunset diagram in the ghost equation given diagrammatically in 
Fig. With the bare four-ghost vertex given in eq. ()97|) and the abbreviations for the 
squared momenta x = p 2 , y\ = (qi) 2 , y 2 = fe) 2 , Z\ = (p — qi) 2 and z 2 = (p — q% — q 2 ) 2 the 
sunset diagram reads 

v- = ffg(i-?) a) 2 U qi f i% (122) 

2(2tt) 8 V2V 2 J ) J xy 1 J y 2 z 2 

The factor 1/x in the first integral stems from the left hand side of the ghost equation. We 
now integrate the inner loop with the help of formula ()108j) and obtain 

U - 512^ \2\ l -2) X ) Y 2 {l-(3)T{2 + 2(3) j dqi l^ (Zl) ' (123) 

where z\ is the total squared momentum flowing through the integrated loop. The second 
integration is done in the same way. We arrive at 

V sun = x w N c9*Z A B* fa / ax ^ r 3(i + /?)r(-3/?-l) 



512vr 4 V2 V 2/7 r 3 (l -/3)r(3 + 3/3) 

:= x 3/3 (f/ , ) sm - (124) 

As each integration step eats up the two squared momenta in the denominators of the integral 
kernels only powers of x to the anomalous dimensions of the dressing functions in the loop 
(here 2>j3 from three ghost propagators) survive. This mechanism works in the same way for 
all diagrams and explains the pattern in the eqs. (|41j) . (J42j) and (J43j) in the main body of the 
paper. 

Next we evaluate the two contributions in the gluon equation needed for the argument 
below eq. ()46|) . The explicit expressions for the kernels of two-loop gluon diagrams are rather 
lengthy but the calculation is done along the same lines as in the ghost sunset diagram above. 
Therefore we just give the final results: 

int 4(T -27 g 4 N 2 c Z, A 4 r(-l-4 ( T)r(l/2- ( r)r(3 ( T)r 2 (l + ( T) v 

V rprprprp OC 



4096 7T 4 T(4 - 3a) T 2 (2 - a) T(3/2 - a) T(4 + 4a) 
2- 4a (-1 + 3a) (10 + a - 66a 2 + 63a 3 ) (5 + 43a + 47a 2 ) 
x 4 - (V') s ^ , (125) 



W sun 7 ^ g'NjZ.C 3 1 T\l + 5)Y{l-38) } 



1536 7T 4 (l + 35)r 3 (2-5)r(3 + 35) 
x 35 (W')ITl- (126) 
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Finally we calculate that part in the dressing diagram of the ghost equation which contains 
the longitudinal part of the gluon propagator for the special case a — 0, 2. These are the 
linear covariant gauges where L(x) — 1 by virtue of the Slavnov- Taylor identity. Replacing 
dressed vertices with bare ones, however, violates this identity We therefore start with the 
general expression, L(x) = Cx p and investigate whether the limit p — > can be performed 
consistently. With the momentum assignments x = p 2 , y = q 2 and z = k 2 = (p — q) 2 the 
longitudinal part of the diagram is given by 



U 



dress _ 

(•_>-)' 

2 



J z 1 xy 



N c g 2 Z 1 \BC f 4 f q^ p^p u q^ 1 p^ q^ 

q l^„l-/3*2-p t* ,,1-/3*2-0 o,-0 v 2-p + 



(2tt) 4 J * y» y i-P z 2- P x y l ~ p z 2 - p y-^z^-P x y-Pz 2 -? 

(128) 

where again the extra factor 1/x stems from the left hand side of the ghost DSE. At this 
stage of the calculation it is not clear whether there are infrared singularities in the limit 
p — > 0. We employ the tensor integrals given in appendix C, use xF(x) = T(l + x) and 
obtain 

dress = N c g 2 Z x \BC 0+p T(2 + /3)T(p)T(-l3-p) p 2 + p/2 

L 16n 2 T {-f3)T{2-p)T{2 + p + {3)f3{2 + f3 + p)' 1 ' 



In the limit p — > this expression is infrared-finite, as lim p ^ T(p)p = 1. We then obtain 

2 ZjXBC 1 
~16^ X 2/3(2 + J) 



ut ress = _ Nc9 ZyXBC j 1 (13Q) 
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